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Abstract. We compute the 0{1/N) correction to the stability critical exponent, w, in the 
Landau-Ginzburg- Wilson model with 0{N) x 0{m) symmetry at the stable chiral fixed point 
and the stable direction at the unstable antichiral fixed point. Several constraints on the 0{1/N) 
coefficients of the four loop perturbative /3-functions are computed. 
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1 Introduction. 



The renormalization group analysis of spin models has proven to be an important tool in un- 
derstanding and predicting critical properties of phase transitions in a variety of phenomena. 
For a recent review see, for example, For instance, the Heisenberg model based on an 0(3) 
nonlinear sigma model or 0(3) (j)^ theory has been widely used to understand ferromagnetic 
phase transitions in nature. Indeed the perturbative field theory techniques in this instance 
have been developed to five loops in MS in standard d-dimensional regularization in cj)'^ theory, 
1^], and equally impressively to six loops in the three dimensional fixed dimension renormaliza- 
tion, 1^, ^. Given such success these methods have been applied to similar models of other 
critical phenomena. Over a period of years various calculations have been carried out in the 
Landau-Ginzburg- Wilson model which essentially is an extension of cp'^ theory where the 0{N) 
symmetry is replaced by an 0{N) x 0(2) symmetry, 0, |9|, 10, 11|. The critical structure in 



this generalised model is richer than the usual 0^ theory in that theoretically there are several 
fixed points over and above the Heisenberg one, depending on the value of N, which are either 
stable or unstable, |1C, However, the properties of the phase transitions in this Landau- 
Ginzburg- Wilson class of models is controversial, g, |ll|, H, H, II3> ID' H- 
First, the two dimensional nonlinear sigma model with the same symmetry is believed to reside 
in the same universality class, ||lO| , Therefore, it ought to be possible to use either model to 
compute useful information on the critical properties of the physically interesting stable chiral 
fixed point. However, it has been pointed out in ||l^ that the results from a d = 2 + e dimen- 
sional study do not match those deduced from the higher dimensional theory. Second, in field 
theoretical calculations such as perturbation theory the phase transitions are regarded as second 
order whilst numerical or Monte Carlo simulations would appear to indicate transitions are first 
order in nature, |^, 16, 17, |l^. Moreover, the particular behaviour depends on the 



value of N though the precise range where this occurs is still undetermined. Further, one recent 
work has suggested an interesting point of view for the origin of the disagreements for N = 2 and 



3 in both two and three dimensions. In [22] it is argued that it is due to the fact that one fiows 
to the stable chiral fixed point along a spiral- like trajectory in contrast to the usual descent. To 
endeavour to clarify the issue the perturbative analysis of the Landau-Ginzburg- Wilson critical 
behaviour has recently been extended to a higher loop order in Previous one and two loop 
computations were carried out in |^ 0, §]. The new three loop MS calculations of the renor- 
malization group functions such as the /3-functions of both coupling constants and anomalous 



dimensions, |11|, have provided more accurate information on the fixed point locations and the 
range of parameters for which they exist and are stable or not. Indeed in this respect the models 
with the more general symmetry group of 0{N) x 0{m) were studied with m only set to m = 2 
at the end, |l^]. Such three loop calculations represent the current perturbative status of the 
model. 

However, it is in principle possible to extend this three loop MS dimensionally regularized 
calculation to the next order, though it will involve a huge number of Feynman diagrams. In 
previous work in the simpler 0{N) models the perturbative computations were complemented 
with large calculations of the same renormalization group functions to several orders in 
powers of 1/A^. For instance, various critical exponents are available at both 0{1/N'^) and 
0{1/N^), p3|, p4|, 25, 26|, as functions of d, with 2 < d < 4, which correspond through the 



critical renormalization group equation with the renormalization group functions. More correctly 
these critical exponents were computed in d-dimensions at the non-trivial Wilson-Fisher fixed 
point of the d-dimensional /^-function which corresponds to the Heisenberg fixed point in three 
dimensions in 0{N) (j)^ theory or the 0{N) nonlinear sigma model which does lie in the same 
universality class. The coefficients of the powers of e in the e-expansion of such exponents 
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in d = 4 — 2e dimensions are in exact agreement with the perturbative coefficients in the 
renormalization group function to the perturbative order they are known, at a particular order 
in More significantly the large critical exponents contain new higher order information 

in the uncomputed coefficients which would therefore assist future perturbative calculations. 
Indeed in the Landau-Ginzburg- Wilson context various critical exponents have already been 
computed in the model with 0{N) x 0{m) symmetry at the two non-trivial fixed points which 
exist in addition to the Heisenberg fixed point, [P, 11]. These are known as the chiral stable. 



(CS), and antichiral unstable, (AU), fixed points. Moreover, the results for the critical exponents 
r/ and v at 0{\/N'^) at both CS and AU are in agreement with the new perturbative results 
of [|ll|. However, to extract the new information encoded in the exponents at four and higher 
loops in these d-dimensional functions in relation to the four dimensional theory one requires 
knowledge of the location of the fixed points to the same loop and large N order. From the 
critical renormalization group equation such information is encoded in the critical exponent 
LO which relates to the critical /3-function slope of the model in the universality class which 
is renormalizable in four dimensions. In the Landau-Ginzburg- Wilson model this has not yet 
been computed at 0{\/N) at either CS or AU. Therefore, it is the purpose of this article to 
rectify this gap and thereby unlock the door to higher order information on the structure of 
the perturbative renormalization group equations such as the /3-function. In 0{N) (j)'^ theory 
this problem has already been resolved at 0(1/A^^), ||2^, where the 0(1/A^) value for oj is 



relatively trivial to establish, [27|, with the elegant machinery of the large N critical point 
method of [psj p^ . However, for the CS and AU Landau-Ginzburg- Wilson fixed points the 
leading order, 0{1/N), analysis is much more involved since one is studying a model with 
two independent coupling constants. Therefore, whilst our calculation also opens the road to 
an 0{1/N'^) computation, it represents a non-trivial example of how one treats the large 
formalism for cu exponents explicitly in a quantum field theory with more than one coupling 
constant which deserves detailed treatment. 

The paper is organised as follows. In section two we recall the background details of the 
model we are interested in and derive explicit expressions for the location of the various fixed 
points from the explicit three loop perturbative results at 0(1/A^) as well as the perturbative 
values of the eigenexponents of the stability matrix at criticality. These are related to the 
exponents uj which we are interested in. Section three is devoted to the development of the large 
N formalism for computing these various to and the explicit d-dimensional expressions are given 
at 0(1/A^). Various concluding remarks are given in section four. 



2 Background. 



The lagrangian of the massless Landau-Ginzburg- Wilson model involves a scalar field with two 
quartic self interaction terms with an 0{N) x 0{m) symmetry and is given by 



L 



ai\2 



+ 



V 

+ 4! 



(2.11 



where l<i<A^, l<a<m, and u and v are the bare coupling constants. As in [^] we rewrite 
(|2.1| ) in order to perform the large N expansion. This involves introducing two auxiliary scalar 
fields one of which, T"^, is a symmetric traceless tensor under 0{m). Thus (2T) is equivalent to 
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where w = u + (m — l)v/m in our notation. If one uses the equations of motion for a and T""^ 
then (|2.1|) is recovered. The coupling constants are defined in the kinetic terms to ensure that 
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the vertices are in the right form for applying the uniqueness technique to compute the large 
Feynman diagrams |2^, . One can understand the fixed point structure of (^]^) by considering 
the /?- functions for each coupling which have been computed to several orders, 11|. At three 
loops these are 
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and 
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1 



(2.4) 



where we have rescaled the coupling constants by a numerical factor to ensure the expressions 
are in the correct format for comparing with the Heisenberg large value for to and the ones 
we compute here. Also we have included the d-dependent terms as we are interested in the fixed 
point structure in d-dimensions. To assist with determining new information on the four loop 
structure of both /3-functions at 0(1/A^) we have introduced parameters, {aj} and for the 
coefficients of the possible terms. 



4 




Figure 1: Renormalization group flow on the {u,v) plane illustrating the Gaussian (G), Heisen- 
berg (H), stable chiral (CS) and unstable antichiral {AU) fixed points. 



By examining the solutions to (3u = and /3„ = 0, ||8|, 11|, several fixed points emerge. First, 
there are the two obvious ones of the Gaussian fixed point, Uc = 0, Vc = 0, and the Heisenberg 
fixed point, Uc 0, Vc = 0. For the latter point setting v = and m = 1 in ( |2.3| ) one recovers 
the usual 0{N) symmetric 0^ theory whose /3-function is known at five loops in MS in four 
dimensions, [0. Indeed in our choice of parametrization of (3.u{u,v) at four loops we used this 
fact to restrict the function to be proportional to v. These two fixed points clearly lie on the 
axes of the {u, v) coupling plane. However, for a range of values of and m there are two 
other fixed points which both have Uc ^ and Vc 7^ 0. One is known as CS and the other AU. 
The ultraviolet renormalization group flow of the four fixed points in the {u, v) plane is shown 
graphically in Figure 1. The range of values for N and m for which such a renormalization 
group flow is present has been detailed in ||8|, for example. However, for the purposes of the 
large N calculation we will require the values of Uc and Vc to several order in e and powers of 
where we take the convention d = 4 — 2e. In |^, 11| the explicit functions of and m were 
presented at two loops. The full expressions at three loops can also be derived but are large, 
|11|, and the full form is not necessary for our purposes. Indeed if we write 
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(2.5) 



for the expansion of the critical couplings in powers of e the four fixed points are determined 
to O(e^) and 0{1/N'^) as follows. For the Gaussian fixed point Uri = Vri = 0. At the Heisenberg 
fixed point un = Q/m, u^i = for r 7^ 1, U12 = — 48/m^, U22 = 108/m^, = — 99/m^, 
U42 = — 7776(ai — a2 — 03 — 04 — «5 — CLQ)/m^ and Vri = 0. At the stable fixed point, CS, 



Mil = 6 , U21 = , U31 

U22 = 18m + 90 , U32 = 

vii = 6 , V21 = , V31 = 

V32 = - 5 (39m + 99) , Vi2 = 



, U41 = 
i(39m + 159) 



U12 = — 6m — 42 
U42 = — 7776ai 

V41 = 0vi2 = — 6m — 24 , V22 = 18m + 54 
777662 . (2.6) 
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Finally, at the unstable fixed point, AU, 
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(2.7) 



These agree to two loops with the expressions given in |8|, |Tl|. As we will be computing the 
critical exponents to which relate to the critical slope of the /3-functions in large we can use 
these values to determine the 0{1/N) form of the critical exponents. As we are working with 
a two coupling model the stability exponents of each fixed point are related to the eigenvalues, 
A/, of the matrix of derivatives, 0,{u,v), evaluated at the appropriate fixed point where 



n(n, v) 



d/3uiu,v) d/3u{u,v) 
du dv 

du dv 



(2. 



For the Gaussian case this is trivial and will not concern us here. For the Heisenberg fixed point 
^l{u,v) becomes triangular because [5v{u,v) has no terms involving only u at any order which 
implies 



du 



Heis 



. 



(2S 



Thus, the critical point eigenvalues in the Heisenberg case ar^ 
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A 
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e + 



mA^ 



12e-10e2 - 13e^ + O(e^) 



+ O 



1 

iV2 



(2.10) 



One of these corresponds to the stable direction in the renormalization group flow, as indicated 
in Figure 1, whilst the other corresponds to the unstable direction which will not concern us here. 
It is for the former for which the corrections at 0(1/A^^) have been computed in d-dimensions. 



*It is worth stressing that our convention, d = 4 — 2e, and the form we take for the /3-functions, (2.3) and 



(2.4), impUes that for a stable direction the eigenexponent is (4 — d)/2 + 0{1/N). This differs from the standard 
result for the stability eigenexponent of (4 — d) + 0{1/N) by a factor of 2 which ought to be taken into account 
when comparing with other calculations. (See, for example, M.) 
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|26|, and we will record the 0{1/N) value later. For the remaining two fixed points the critical 
matrix Q.{uc,Vc) is not triangular. However, computing the eigenvalues at criticality of ( |2.8| ) we 
find 

3(m2 + 4m + 7)e^ (13m^ + 72m^ + 202m'^ + 216m + 25) 3 
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(2.11) 



where the sign of the 0{e) terms relates to the stability property of the fixed point when viewed 
from the ultraviolet renormalization group flow of Figure 1. For the exponents corresponding to 
the stable direction we have included the 0{e^) terms which depend on the unknown parameters 



of the 0{1/N) four loop /3-functions of ( |2.3D and (2.4) and which will be constrained by our 
0{1/N) critical exponents. 



3 Large N formalism. 

To compute the same critical eigenexponents in the large formalism in d-dimensions one 
follows the programme of p3| , p4| where the appropriate Schwinger Dyson equations are analysed 
at the respective d-dimensional fixed points of the theory. At this point the propagators scale 
asymptotically to power law behaviour. In other words, in coordinate space, we have 

^a^,bJ^^^ = , a{x) ~ ^ , a^'^'^x) = X-'''^aT{x) (3.1) 



where 



A C 
'^(^) ~ ' ''Tix) ~ (3.2) 
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and 



1 / o 

2 \ m 



This group theory factor satisfies the projector property 

j^ab,pq j^pq,cd j^ab,cd 



(3.3) 



(3.4) 



which imphes that in labehng the internal indices on a diagram one only needs to put indices 
on the (/(-field lines. The quantities A, B and C are the x-independent amplitudes of the fields 
and the exponents, a, /3 and 7 in our notation, are related to the wave function anomalous 
dimension r/ by 
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1 + ir/ 
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(3.5) 



where x ^^'^ XT are the respective anomalous dimensions of the vertices involving a and T'^'' 
and d = In |ll| r] was computed at 0(1/ N'^) at both fixed points CS and AU. For the stable 
chiral one both the a and T"^ fields propagate and couple. However, at the Heisenberg fixed 
point only the a field propagates since clearly the T"* field is absent in the usual formulation of 
0^ theory. Interestingly at the unstable antichiral fixed point the opposite situation emerges in 
that only the T"^ field is present and the a field is omitted from the calculation of the large 
exponents, [0]. Since the amplitudes appear in the calculations in the combinations z = A^B 
and y = A^C throughout and as they will be required for computing our uj exponents we have 
determined their values at 0(1/A^). For reference, at CS they are 



2T{2n - 2) 



mr(2 - ^)r(^ - 2) 
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and at AU we have 
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Moreover, for completeness the exponent r/ at 0{1/N) is given by 
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Any subsequent exponent at either fixed point will be expressed in terms of their respective 
value for r]i. 

Since the exponents loj relate to corrections to scaling then to compute them one considers 
corrections to the asymptotic scaling forms (p^), [ 
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24, 26 1 . In coordinate space we take 
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I + B'(x 



1 + C'{x 



(3.10) 
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where A' , B' and C are the x-independent correction to scahng amphtudes whose values are not 
important here. In addition to ( 3.10I ) one requires the scahng form of the inverse propagators 
which are determined by inverting the Fourier transform of ( p. 10 ) in momentum space. Thus 



a-\x) 



p{a) 



(x2)2m-/3S 



(x2)2m-7C 

where the functions p{x) and q{x, y) are defined by 
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(3.11) 



(3.12) 



with a{x) = T{fi — x)/T{x). Further, in our notation each exponent ujj in the stable direction 
will have the 1/A^ expansion 

oo 



to 



(/^-2) + E 



(3.13) 



i=l 



and are therefore related to the eigenexponents of ^}{uc,Vc) by a;/ = — A/ respectively in our 
conventions. Moreover, as noted before our choice of uj differs from the usual definition by a 
factor of (-1/2). 



= f ^ 




Figure 2: Schwinger Dyson equation for cp field at 0(1 /N). 

To solve for uJi one substitutes the asymptotic scaling values into the lines of the Schwinger 
Dyson equations of each of the fields. These are illustrated in Figures 2 and 3. In the latter 
we have only indicated the possible topologies that arise due to the nature of the expansion in 
powers of 1/A^ and each wiggly line represents all possible combinations of a and T"* fields. If 
we ignore for the moment the two and three loop corrections in Figure 3 we find the equations 
are represented by 
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+ —Nmz 
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(7,cj)(x2)'^C7' + ^Ny 1 + 2{x^)'^A 



1 + 2(x'^rA' 



„2\LU A' 



(3.14) 



at the CS fixed point which we consider first for illustration. These equations decouple on 
dimensional grounds into a set which determine r/ at CS and a set involving lu. The former lead 
to the values for rji, z\ and y\ quoted above. For the latter set of equations for consistency 
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the determinant of the 3x3 matrix defined with respect to the basis vector {A', B' , C'} has to 
vanish which naively leads to the equation 







(m - l)(m + 2] 



m 



yqi(3,uj) - q{-f,uj) 



2z 



(l + g(a,u;))g(/5,cu) 



z + 



(m - l)(m + 2) 



2m 



(3.15) 

This is similar to the consistency equation which determines the exponent denoted by A in |p3| , p^ ] 
where the correction to scaling exponent has canonical dimension Aq = /x — 1. However, for 
this value q{a,Xo) with a = (/x — 1) or 2 canonically is always an 0(1) quantity. By contrast, 
for the four dimensional model q{a,uJo) with ujq = {fi — 2) and a = (// — 1) or 2 it has values 
which are 0{N) and 0{1/N) respectively and therefore affects the consistency equation. For 
instance, whilst the product q{a, uj)q{f3 , uj)q{'y , ca) in ( |3l5| ) will be 0(l/iV), this equation, 
would give an incorrect value for the exponents since the contributions in ( |3.15| ) from and 
are of the same order in 1/A^ as those from the two and three loop correction graphs of 
Figure 3 which we had naively omitted. Thus to determine wi correctly one must include the 
contributions from these diagrams to the B' and C parts of the Schwinger Dyson representation 
in the asymptotic approach to the fixed points. Hence, the a and T"'' equations in ( 3.14| ) are 
modified to 
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(3.16) 



In these equations to ensure the correct contribution to the uj Schwinger Dyson equation after 
decoupling one correction term, (x^)'^, is on one internal a or line. There are corrections to 
the ^ lines but these only contribute to uj2 after examining where they appear in the consistency 
determinant. This feature of having to include higher order graphs to obtain the correct bj\ is 
not peculiar to the CS fixed point as it already occurs at the Heisenberg point, [^]. With these 
additional diagrams the correct consistency equation is given by setting the determinant of the 
matrix 



(l + (7(a,a;)) 



2 + 
2 
2 



(n^_l)(n^+2) - 
2m y 



q(/5, u)) + zHi + 2mNz'^U2 
z{Ui + 2A^yn2) 



(m-l)(m+2) 
2m y 

n(in±^y^m + 2NyU2) 



2m 

q{j,uj) + 

, (m-2) 



2m 



yUi + 2NyzIi2 



(3.17) 



to zero. 



Therefore, to solve for the respective w's from ( |3.17 ) the values of these additional diagrams 
must be computed. As the graphs are similar to those used in the 0{N) cfA calculation of oj if 
the group theory of each graph is suppressed we merely quote the d-dimensional values for the 
respective two and three loop graphs of Figure 3. They are, ignoring symmetry factors. 



Hi = [v{2,^l-l,^^-l)f 



Ho 



[i/(2,/i-l,/x-l)]2zy(l,2,2^-3)z^(4-^,/x-l,2/u-3) (3.18) 
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Figure 3: Basic topologies for the corrections to the a and T"^ Schwinger Dyson equations to 
determine a; at 0{1/N). 



where iy{x,y,z) = a{x)a{y)a{z). These were calculated using the method of uniqueness of 



which was extended from three dimensions to d-dimensions in |23, 24 1. Futher, since the 0{1/N) 
values of q{(3,uj) and q{'~f,uj) now need to be included it transpires that the 0{1/N) expression 
for the vertex anomalous dimensions, x xt, are required. This is due to 



/(Aw) 



a(4 - ^)r(^) [uji -r]i- xi] 
a{2)a{2-fi) N 



+ O 



1 

W 



(3.19) 



We have computed each vertex anomalous dimension at each fixed point and find 



Xt,i 



^AU 



CS 



^(4/i - 5)?7^ 

(/x-2) 
/i[(2/i - 3)m + (4/i 



CS 



(^-2)(m + l) 
-2)[(m + 4)(2^ 



3) + lU 



AU 



[m 



l)(m + 2)(/i-2) 



(3.20) 



The expression for x^^ is formally the same as that for the 0{1/N) vertex dimension at the 
Heisenberg fixed point. For each exponent we have computed the value by applying the tech- 
nique of [^] of large N critical point renormalization of 3-point functions to each vertex at the 
appropriate fixed point. For x we have checked that the value agrees with that given by the 
scaling law which emerges from the theory which is in the same universality class as (|2.2| ). This 
is believed to be the 0{N) x 0{m) two dimensional nonlinear sigma model where instead of the 
term quadratic in a of ( |2.2D one has a term linear in a where its coupling constant is related to 
the critical exponent v. However, in such a model it is clear from group theory that there can 
be no linear term in T""^ and therefore both expressions for xt cannot be derived from a scaling 
law but only direct calculation. 

With these values we can now determine the solution for the consistency equation for each 
fixed point. For CS since the matrix is 3 x 3 two values for u emerge. These are 



0^2^ = (/i-2) + 



(2/i - 1)??P^ 
2(m + 1)(^ - 2)N 



m{fi - 1)(^ - 4) + (2/^2 -7^ + 4) 



(m^ - 1)(^ - 1)2 + 2(m - l)(2/i - 3)(/i - 1) + (5^ 



+ O 



1 

(3.21 



where the ± subscript refers to the sign in front of the discriminant. Both are perfectly acceptable 
solutions since one corresponds to one stable direction at CS and the other to the eigenexponent 
from the second direction. 

For AU the derivation of the consistency equation follows the same pattern as that for CS 
in that two and three loop graphs have also to be included due to the large N counting. The 
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consistency equation itself can be deduced from ( 3.17 ) by deleting the second row and column 
from the matrix and setting z = in the remaining entries. Since only the T"-^ field is present 
only one value for io emerges which is 



AU 



- 3/i - 1 + 



H{m - 2)[2/i - 5 



(m - l)(m + 2) 



2(m + 4)(2^ - 3)]] 77j^U 



+ O 



1 

W 

(3.22) 

which corresponds to the eigenexponent of the stable direction at AU. For completeness we note 
that at the Heisenberg fixed point for the stable direction we have, |27, p^ ], 

4(2// - l)2r(2/i - 2) 



CO 



Heis 



(/i-2) 



r(2 - i2)r{n - i)r{fi - 2)r{n + i)mN 



(3.23) 



which is deduced by deleting the third row and column of (|3.17| ) and setting y = 0. 

In order to check the correctness of each result one can set d = 4 — 2e, or /i = 2 — e, and 
expand each exponent to O(e^) and compare with the values of explicit perturbation theory, 
( 2.11 ). We find total agreement for all three exponents and therefore regard ( 3.21 ) and ( |3.22 ) 



as correct. Indeed equipped with ( p.21|) and ( 3.22 ) we can now reverse the check argument and 
use the information contained in the d-dimensional exponents to determine constraints on the 
structure of the four loop MS corrections to ( |2.3D and (|2.4| ). First, we record the e expansion 
of each of the exponents in the stable directions which are consistent with ( 2.11| ) . From ( |3.21| ) , 
(|0^ ) and (^) we have 



cs 



CS 

-.1 



AU 
+,1 



Heis 
+,1 



3(m2 + 4m + 7) 



(13m'^ + 72m3 + 202m'^ + 216m + 25) 



(m + 1) 



2(m + 1)3 



+ 



(3m6 + lOm^ + 21m^ - 36m^ - 387m^ - 326m + 395) 



4(m + 1)5 



+ 0{e 



6(m + l)e 



(13m2 + 26m + 25) ^ (3m^ + 12m3 + 82m2 + 116m - 5) 



(m + 1) 



+ 



2(m + 1)3 



_^ [48(m + 1)^C(3) + 3m6 + ISm^ + 29m'' + 76m3 + 397m2 + 322m - 397] ^4 



+ 0(e^ 



1 

m 
1 

m 



(3m2 + 9m - 34) 



4(m + 1)5 



-(13m2 + 33m- 162)e3 + ^(3m2 - 5m - 54)e^ 
2 4 



18e^ 



33e^ 



(3.24) 



Next we compare the O(e^) coefficients with those of the explicit perturbative expansion which 
have been parametrized by {aj} and {bi}. From the Heisenberg exponent we have 



ai 



a2 



as 



04 



05 



06 



5m 
7776 



(3.25) 



which agrees with the explicit four loop (j)^ computation when one sets m = 1 and determines 
the 0{1/N) coefficient of the term of f5u{u,v). From the CS fixed point we have 



2ai + (m - 1)62 



(3m + 7)(m — 3) 
15552 



(3.26) 
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and 



(m + l)[a2 + 2a3 + 804 + 4a5 + Soe - 63 - 264 - 865 - 466] 
+ (2m — l)(m — l)ai — 2m{m — 2)62 
C(3) ^ (3m2 + 6m + 19) 



108(m + l)3 5184(m + l 

The final constraint arises from the stable direction at the AU fixed point which gives 

(m - 1)^62 + (m - 1)^63 + (2m - l)(m - 1)^64 + {3m^ - 8m + l)(m - 1)65 
+ (2m - l)(2m2 - 2m + 1)65 
[8m^ - 5m - 54] 
15552 



(3.28) 



One can now deduce the values for the exponents in the four stable directions in three 
dimensions. These are, with our conventions, 



CS _ 1 4(m + 4) /I 

- " 2 + ^ 
CS _ 1 16(m + l) /I 



AU _ _ 1 4K+4m-8] ^ 
Heis 1 32 „ / 1 



For CS the corrections both have the same sign and interestingly neither involves a square root 
which appears in the d-dimensional expression. For the specific case of m = 2 we have from 



,,AU 



, ,Heis 
OJ, 



- -- of — 

1 16 ^fl 



m=2 2 7r2iV KN"^ 

m=2 ~ 2 37r2iV ViV2 



+ ^ + O • (3-30) 



m=2 2 37r2iV ViV2 



With these values we can comment on the possible breakdown of stabilty in this model. In 
our computations so far have relied on the fact that the stability picture for the fixed points 
represented in Figure 1 is valid for the full range of in the large N expansion. However, 
it has been suggested that for certain values of N this scenario may be different and that the 
underlying assumption of the existence of a second order phase transition, around which the 
large N critical point formalism is built, could break down. Indeed there is some controversy in 
this model in three dimensions about the existence of CS and whether the phase transition is 
first or second order for relatively low values of N . Moreover, the precise value of N where the 
order changes has not been determined consistently from different methods. (A recent review 
is given in |]l^.) For instance, a value has been found for this critical value of N by using 
d = 4 — 2e dimensional perturbation theory and it was extended to three dimensions using 
standard resummation techniques, pO], giving Nc = 3.39. By contrast, Monte Carlo methods 
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and another e-expansion extraction have suggested Nc <2, 0, ^. With the large corrections 
( 3.30| ) we can naively examine the range of values of N for which either CS stability exponent 
changes sign. For to^^ this will occur when Nc = 3.24 whilst a;9^ changes sign when 



=2 



=2 



Nc = 1.62. The former would suggest a critical in a similar range to that of |30|]. However, 
these remarks ought to be qualified with various observations. First, we have only computed 
the 0{1/N) correction to stability where is assumed to be large. Therefore, one has to ask 
whether the approximation will still be valid for such a low value of A^. Second, the nature of the 
large A^ expansion is a reordering of perturbation theory such that a certain class of diagrams 
are summed first. Therefore, if one could compute to all orders one would reproduce ordinary 
perturbation theory and so obtaining a value for Nc in three dimensions which is not inconsistent 
with the resummed value determined from several loop orders in ordinary perturbation theory 
would seem only to reinforce that particular value. In other words if non-perturbative effects 
become significant at CS for low values of A^ to affect the precise location of Nc these will have 
been omitted in perturbation theory. Third, for our value of Nc we have naively assumed the 
large A^ series is convergent and therefore that our simple assumption that when the 0{1/N) 
correction exceeds 1/2 the character of the fixed point changes is valid. Only a higher order 
calculation would resolve this. 



4 Discussion. 



We have computed the correction to scaling exponents co in all the stable directions of the 
Landau-Ginzburg- Wilson model with 0{N) x 0{m) symmetry at 0(1/A^) in d-dimensions. 
This allows one to extract information in all the available large A^ exponents of |11] at 0(1/A^) 
in relation to the four dimensional theory in the same way that the two dimensional critical 
slope information contained in the exponent i' does for the underlying two dimensional theory. 
It is also worth stressing again, [0], that from the point of view of the large A^ formalism a 
consistent picture emerges in terms of the active fields of the theory formulated in terms of 
the auxiliary fields a and T°'^ of (^.2|). At the Heisenberg and AU fixed points only a and 
rpab respectively propagate which corresponds in the large A^ formalism developed here to one 
stable direction and hence only one eigenexponent emerged. However, at the only fully stable 
fixed point both fields are relevant leading to two independent stability exponents. This is a 
natural way to picture this particular model which we assume persists to higher orders in large 
A^. Moreover, our consistency with three loop MS perturbation theory provides an important 
internal cross check on the values of quantities, such as the vertex anomalous dimensions, which 
we had to compute en route to our expressions for Wj. Further, the information contained within 
the expressions ( [3.21 ), ( |3.22| ) and ( [3.23 ) will provide important checks on any future explicit 
four loop MS perturbative calculations which would improve the accuracy of the numerical 



estimates deduced from (2.3), ( |2.4| ) and other renormalization group functions. Such four loop 
calculations are certainly viable since five loop results are available in MS in ordinary 0{N) 
(j)'^ theory. For example, the integration routines for the four loop Feynman diagrams have 
already been constructed. However, one can also attack this problem from the large A'^ point 
of view. For instance, we have demonstrated the elegance of the formalism to produce the 
critical eigenexponents at 0{1/N). However, this machinery has already been extended in [^] 
to compute UJ2 in the Heisenberg model. Therefore, we would expect there to be no serious 
obstacles to extending the present computation to find cj2'^ and ^2^^. For example, the values 
of the underlying three, four, five and six loop Feynman diagrams which are analogous to the 
corrections of Figure 3 for the 0(1/A^^) calculation have been determined. We hope to return 
to this in a future article. 
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